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Abstract 

The origin of Ultra-High Energy (E Z 10 20 eV) Cosmic Rays (UHECR) remains 
mysterious. I discuss "top-down" models, where UHECR originate from the decay of 
very massive, long-lived particles. I summarize the calculation of the spectrum of decay 
products, discuss possible problems with this scenario, and describe ways to test it by 
searching for very energetic neutrinos and neutralinos. 



The observation [1] of CR events with E > 10 20 EeV poses two problems: 

• The energy problem: Few, if any, known celestial objects possess sufficiently strong 
electromagnetic fields extended over a sufficiently large length to accelerate protons to 
such energies. 

• The propagation problem: Protons with E > 5 ■ 10 19 eV lose energy by photopro- 
ducing pions on CMB photons [2]. This should lead to a sharp drop in the spectrum 
if sources are distributed homogeneously (the "GZK cut-off"). Similarly, heavier nuclei 
lose energy through photo-dissociation, and photons get converted into e + e" pairs [3]. 

The second problem aggravates the first one: The GZK effect implies that sources of UHECR 
should be within ~ 50 Mpc. Even charged particles (in particular, protons) should travel 
along rather straight lines over such distances, once E > 10 19 eV. However, there are no 
known nearby sources in the UHECR arrival directions. 

In top-down models [4] the energy problem is solved by postulating the existence of "X" 
particles with mass mx > 10 12 GeV and lifetime Tx > 10 10 yrs ~ t\j. The decay of such 
objects clearly can generate particles of the required energy. Moreover, several mechanisms for 
producing a sufficient abundance of X particles in the very early Universe have been proposed 
[5] . Finally, X particles can be very long-lived either by embedding them in topological defects 
[4, 3] , or by greatly suppressing their couplings to lighter particles [6] . Essentially X particles 
are batteries storing energy from an earlier, much more violent epoch until the present time. 
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The propagation problem can also be solved in this scenario. If X particles are distributed 
more or less evenly throughout the Universe, as expected e.g. in topological defect models [3], 
a spectral break at the GZK energy is expected, in agreement with the HiRes (but not the 
AGASA) data, but there might be some "nearby" sources not associated with known celestial 
objects. Recall that HiRes does see some post-GZK events, which implies the existence of 
"nearby" sources (unless Lorentz invariance is violated [7].) Alternatively, if X particles move 
freely under the influence of gravity, an enhancement in their density by a factor ~ 10 5 is 
expected in the halo of our own galaxy. In this case almost all observed UHECR events would 
be of "local" origin, and no sharp drop-off would be expected at the GZK energy, in agreement 
with the AGASA (but not the HiRes) data [1]. In both cases the spectrum should cut off at 
E = mx/2. 

In order to test this idea quantitatively, one first has to compute the spectrum of X decay 
products. The very existence of an energy scale m x strongly indicates the existence of weak- 
scale supersymmetry (SUSY), in order to stabilize the hierarchy between mx and the weak 
scale (~ 100 GeV) against radiative corrections [8]. One then expects primary X decays to 
produce comparable numbers of ordinary particles and their superpartners. This primary X 
decay will start a parton cascade, similar to Z — > qq decays at LEP starting a QCD cascade 
leading to a multi-parton state. There are a couple of important differences, though, related 
to the fact that the particles produced in primary X decay will have time-like virtuality near 
mx- First, at such scales all three gauge interactions have comparable strengths, and should 
therefore be treated on equal footing. Secondly, at such scales superparticles can be treated 
as massless, and can thus be produced during the cascade, along with quarks, gluons, leptons 
etc. 

Once the cascade virtuality scale drops below the mass scale of SUSY particles m S usY < 1 
TeV, superparticles (as well as top quarks and electroweak gauge and Higgs bosons) decouple 
from the cascade and decay (often starting their own little QCD cascades). At virtuality 
< mgusY only strong interactions need to be kept in the description. Finally the virtuality 
scale reaches Qhad — 1 GeV, where partons hadronize into baryons and mesons, many of which 
in turn decay (along with the heavy leptons). 

This cascade started by X decays has been treated in [9] using the language of fragmenta- 
tion functions, which allow relatively easy calculation of the single-particle inclusive spectra^, 
exactly the quantities needed in our case (since at most one particle produced in a given X 
decay will be detected on Earth). Not surprisingly, the high-energy end of the resulting spec- 
trum depends strongly on the most important primary X decay mode(s). If X decays into 
ordinary particles, most of the energy will be carried away by neutrinos and photons; protons 
carry at most 10%. If the primary decay products are superparticles, LSPs often carry the 
largest fraction of the released energy; hundreds of LSPs are produced even in X — > qq de- 
cays. On the other hand, the shape of the proton, photon, electron, and electron and muon 
(anti) neutrino spectra at energies < O.Olmx are essentially independent of the primary X 
modes, and depend only logarithmically on mx if they are expressed in terms of the scaling 
variable x = 2£ partic i c /m x . 

Top-down models face several challenges: 

• For all primary X decay modes that have been investigated [including the whole cata- 
logue of (s)particles that exist in the MSSM] the 7 flux is significantly higher than the 
proton flux at source. This is problematic, since various independent investigations of 

^Our code computing the X decay spectrum is publicly available [10] 
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Haverah Park, Fly's Eye and AGASA data indicate that relatively few, if any, UHECR 
events are due to photons [11]. If most sources are at cosmological distances, propagation 
effects increase the p/7 ratio. For sources in the halo of our galaxy propagation should 
be significant only if the galactic radio background has been underestimated by at least 
a factor of 10. Moreover, absorbing UHE photons has its own problem: It initiates an 
electromagnetic cascade, which eventually produces photons in the 100 MeV to 1 TeV 
energy range probed by the EGRET satellite. The upper bound on the isotropic flux of 
such photons imposes strong constraints on top-down models [12]. This is probably the 
most serious current problem for top-down models. Possible excuses are the rather poor 
statistics at very high energy, the fact that constraining the observed p/7 ratio heavily 
relies on MC studies, or the fact that the crucial EGRET bound has not been checked 
independently. 

• Precision measurements of the CMB anisotropy [13] indicate that topological defects 
play little or no role in structure formation in the Universe. It is not easy to square this 
with the requirement that the (symmetry breaking) scale associated with the defects 
has to exceed 10 12 GeV in order to explain the UHECR. Freely moving particles thus 
seem the more likely option. As noted earlier, they should be concentrated in the halo 
of our galaxy. Since we are not in the center of our galaxy, one expects an anisotropy 
in the UHECR arrival directions in this "decaying dark matter" (DDM) model [14], in 
particular a peak roughly in the direction of the galactic center, where the DM density 
should be highest. This location is not visible from the Northern hemisphere, but could 
be observed by the Australian SUGAR array. The data from this experiment show no 
such peak [15]. However, the size and exact location of the expected peak depends on 
details of the DM distribution near the galactic center, which is poorly understood, and 
also on the extension of the dark halo of our galaxy, which is unknown. Even taking 
standard halo models at face value, SUGRA data are compatible at the 10% level with 
all events above 6 • 10 19 eV being due to DDM. On the other hand, if the Auger array 
in Argentina also fails to see any enhancement towards the galactic center, the DDM 
variety of top-down models would be in serious difficulties. 

• The AGASA data above 10 19 eV show quite pronounced (at the 4 to 5a level) clustering 
on small scales [16]. The HiRes collaboration is currently not claiming such clustering, 
but also does not seem to be willing to state (yet another) discrepancy with AGASA. 
Although some small-scale anisotropy can be expected in DDM models with clumpy 
halo, the dumpiness required to explain the AGASA data does not seem to be compatible 
with current halo models [17]. Note that here the structure of the halo in our "vicinity" 
is most important, which is supposed to be fairly well understood. However, even in the 
AGASA data the evidence for clustering becomes weak if one restricts oneself to events 
(well) above the GZK cut-off; this is opposite to what one expects in bottom-up models, 
where the number of relevant sources should decrease rapidly beyond the GZK cut-off. 
In any case, the present data on small-scale clustering do not exclude all UHECR with 
E ^ 10 20 eV being of top-down origin. Similar remarks apply for the claimed correlation 
of UHECR with E ~(a few) times 10 19 eV with certain BL Lac's [18]. 

• The last two points indicate that many, perhaps even most, events right at and below 
the GZK cut-off should not be of top-down origin, whereas the energy problem indicates 
that this should be the dominant source at E > 10 20 eV. Shouldn't one expect some sort 
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of spectral feature in the transition region? Generally speaking, Yes, but this argument 
is difficult to make more quantitative as long as the origin of CR's with E ~ 10 19 eV is 
not understood in detail. Note also that the AGASA data are compatible with a spectral 
break (a hardening of the spectrum) in this energy range. 

It should be noted that top-down models are more vulnerable than "generic" bottom-up 
models because they make detailed predictions for the spectrum at source (for given X decay 
mode), and in case of DDM models also about the distribution of sources. In contrast, bottom- 
up "models" often simply assume that protons are injected with a given spectrum (usually 
a power law with a cut-off), without bothering to explain how the protons were accelerated. 
Moreover, often a continuous distribution of sources is assumed (e.g. in the "model" adapted 
by HiRes to explain their data [1]), even though the number of possible bottom-up sources 
within a few GZK interaction lengths is certain to be quite small. To my mind this is only a 
fancy parameterization of the data, not a physical model. 

As discussed above, the more predictive, and currently favored, DDM variety of top-down 
models makes predictions that can be tested by Auger. Other tests of top-down models are 
based on predictions for neutrino[19] and neutralino (LSP) fluxes. These signals would be even 
bigger (by a factor ~ 10) if most sources were at cosmological distances, since one would then 
require a larger flux at source to explain the UHECR data. In DDM scenarios, normalizing 
the observed UHECR flux to the proton flux only, one expects the following signals: * 

• Neutrinos [20]: 1-30 (mostly v^) events/yr above 100 TeV in IceCube, 0.4-4 (mostly v e ) 
events/yr in RICE, and 0.3-3 (mostly v T ) events/yr in Auger. The ranges correspond to 
variations in rrtx and the primary X decay mode, with smaller mx and less hadronic X 
decays giving higher rates. However, less hadronic decays also give a larger 7/p ratio at 
source, and hence a bigger problem with either the data [11] favoring protons as UHECR 
primaries, or with the EGRET bound [12]. However, most other UHECR models also 
lead to detectable neutrino fluxes at E > 100 TeV, where the atmospheric background 
becomes negligible [12]. The neutrino flux should thus give valuable quantitative tests, 
but does not distinguish models qualitatively. 

• Neutralinos [21]: Only top-down models can produce a detectable flux of UHE neu- 
tralinos, since only here the flux of particles and sparticles at source are comparable. 
Note that LSPs are expected to have 10-100 times smaller interaction cross section with 
matter than neutrinos [22]. There thus exists a range of energies (between a few 10 8 
GeV and a few 10 10 GeV) where neutralinos can still cross the Earth more or less unmo- 
lested whereas neutrinos get stuck. Looking for UHE LSPs thus boils down to looking 
for clearly upgoing events in this energy window. Unfortunately the expected event rate 
is marginal even for (first-generation) space based experiments like EUSO, but at least 
in principle this signal is detectable. 

In summary, CR events with E > 10 20 eV may well originate from the decay of very 
massive particles; however, this contribution to the CR spectrum should be sub-dominant for 
E < 5 • 10 19 eV. There may therefore exist two distinct UHECR puzzles. The component (s) at 
E < 5 ■ 10 19 eV would have to come from some acceleration mechanism. Crucial tests of this 
picture should come from better measurements of the UHECR spectrum, and from analyses 

*Rates are typically three times smaller if the sum of photon and proton fluxes is normalized to the UHECR 
flux. 
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of the clustering of UHECR, and/or their correlation with known celestial objects (e.g. the 
center of the Milky Way). A complete understanding will probably require studies of the 
UHECR sources using penetrating particles. In this picture, both bottom-up and top-down 
sources are likely to contribute to the spectrum of neutrinos at E u > 100 TeV; a detailed 
measurement of the spectrum of these neutrinos, while difficult, should then be rewarding. 
Detecting UHE neutralinos would constitute a "smoking gun" signature for top-down models. 
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